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Many body effects on conduction and diffusion of electrons and holes in a semiconductor quantum 
well axe studied using a microscopic theory. The roles played by the screened Hartree-Fock (SH ) 
terms and the scattering terms are examined. It is found that the electron and hole conductivities 
depend only on the scattering terms, while the two-component electron-hole diffusion coefficients 
depend on both the SHF part and the scattering part. We show that, in the limit of the ambipo- 
lar diffusion approximation, however, the diffusion coefficients for carrier density and temperature 
are independent of electron-hole scattering. In particular, we found that the SHF terms lead to 
a reduction of density-diffusion coefficients and an increase in temperature-diffusion coefficients. 

Such a reduction or increase is explained in terms of a density- and temperature-dependent energy 
landscape created by the bandgap renormalization. 
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Our understanding of Coulomb interaction in an op- 
tically excited semiconductor has been greatly enriched 
recently thanks to the extensive theoretical 1-3 and ex- 
perimental investigations 4 over the past few decades. As 
far as optical properties in highly excited semiconductors 
(with high carrier density) are concerned, many-body ef- 
fects manifest themselves in two important ways: the 
renormalization of the single-particle energies and the fi- 
nite life time of such renormalized single-particle states. 
Though the quantitative calculation of these two quan- 
tities is still a topic of current research, the qualitative 
difference of these two features seem to be quite clear. For 
the sake of convenience of the presentation in this paper, 
we classify the two types of terms into coherent and in- 
coherent (or scattering) parts, since the renormalization 
of single-particle state changes the resonance frequency, 
while scattering leads to a decay of coherence or oscil- 
lation. For a spatially uniform system, these two types 
of effects manifest themselves in the linear optical spec- 
trum and have received extensive attention in the past 
decades. The manifestation of the many-body effects in 
a spatially non-uniform system has not received compa- 
rable attention. It is especially interesting to see, as we 
will show later, that these two parts of many-body inter- 
action play different roles in the conduction and diffusion 
processes in a spatially non-uniform semiconductor. 

In this communication, we report on our recent theo- 
retical study on the effects of many-body interaction in 
a spatially inhomogeneous system. The starting point of 
our investigation is the set of Boltzmann-Bloch-Poisson 
equations which contains many-body interactions in the 
spirit of Refs. 1-3 and 5. Namely, the coherent part 
is treated within the screened Hartree-Fock (SHF) ap- 
proximation, while the scatterings are treated within the 
second Born approximation. After following the stan- 
dard moment equation approach 6 by assuming the quasi- 
equilibrium distribution of carriers, a set of coupled dif- 
fusion equations for carrier densities and temperatures 


can be derived 7 and given as follows: 
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where N a , W a , u a , and T a stand for density, thermal en- 
ergy, drift velocity, and temperature of electrons (a = e) 
and holes(o = h). j and are transform Jacobians 
which relate thermal energies to temperatures. They are 
defined in Ref. 7. For the purpose of this paper, it suffices 
to say that they are functions of densities and temper- 
atures. In (1) and (2), s represent generation and 
recombination of carriers due to pumping and optical 
transitions, while R % : s represent the corresponding heat 
sources or sinks. The density and thermal currents in 
equations (1) and (2) are defined as in Ref. 7: 
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The density currents can be written in terms of gradients 
of the four macroscopic variables (X = N e ,N h ,T e ,T h ) 
and the electrical potential, as 7 
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where we have introduced various diffusion coefficients 
and conductivities 7 (a//?): 
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FIG. 1: Diffusion coefficient Dnn (b) and its relative change 
SDnn (a) versus carrier density at three temperatures as in- 
dicated. Solid and dashed curves in Fig. 1(b) are for Dnn 
and Dfv/V) respectively. 
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FIG. 2: Diffusion coefficient Dtn (b) and its relative change 
5 Dtn (a) versus carrier density at three temperatures as in- 
dicated. Solid and dashed curves in Fig. 2(b) are for Dtn 
and Dtn j> respectively. 


D X a T a — Pa (1 + 77a)(C*£ a + ^T a ) + j > (&) 

D N c T 0 — Ha j^(l + Tj a )H^> 3 + Hy# + Cy 0 j , (9) 

^ - fi a e 2 [N a (l + T) a )-N 0 ] . (10) 

Following two shorthand notations have been introduced: 


C% = d x W a , (11) 

H% = N a d x Se a . (12) 

Obviously, the 0 %' s are specific heats of a certain kind, 
which represent the contribution from free electron and 
hole gases, while H# 5 s are the contributions due to many- 
body interaction, as they are proportional to the deriva- 
tives of self-energy renormalization 2 (<5e Q ) with respect to 
densities or temperatures. While equations (6)-(9) define 
diffusion coefficients in the density currents (noting the 
first index of the coefficients being N a ), the correspond- 
ing diffusion coefficients in the thermal currents are de- 
fined through the relation between the density currents 
and thermal currents «/£, Eq. (4). Factors p a and 
in Eqs. (6)— (10) are defined as follows: 
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to carrier-LO (c-LO) phonon and due to electron-hole 
(e-h) scattering, respectively. These rates are defined 
microscopically in Ref. 7. 

Several general features can be readily observed from 
the expressions in Eqs. (6)-(10): First, all diffusio n coe f- 
ficients and conductivities depend on momentum relax- 
ation rates due to e-h scattering through factors p a and 
r/ a . On the other hand, the coherent part of the many- 
body interaction (through H £) only enters the diffusion 
coefficients, but not the conductivities. Second, it is in- 
teresting to consider the ambipolar diffusion coefficients, 
as commonly defined by setting equal the density cur- 
rents of electrons and holes. The ambipolar currents are 
now written as 7 
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where a,/? € {e, h} | a ^ (3. Moreover, y£ 0 and y° h are where X = N,T as we also assumed the densities and 

the relaxation rates of the a component momentum due temperatures to be the same for the two components. 
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The identity in equation (17) defines the free-carrier dif- 
fusion coefficient D° NX and the corresponding many-body 
correction A D^x- We observe that the e-h scattering 
rate yf h (the incoherent part of the many-body effects) 
disappears from the ambipolar diffusion coefficients com- 
pletely. Only c-LO phonon scattering rates 7 £ 0 ’s remain. 
This means that the ambipolar diffusion coefficients de- 
pend only on the weighted sum of the scattering rates of 
electrons and holes with LO phonons. The coherent part 
of the many-body interactions, however, remains present 
in the ambipolar diffusion coefficients, as expressed by 
the second term A Dnx and A Djx in (17) and (18), 
respectively. The absence of the e-h scattering and the 
remaining presence of the coherent part of the Coulomb 
interaction in the ambipolar diffusion coefficients clearly 
illustrate different roles played by the two aspects of the 
same Coulomb interaction. We note that the absence 
of the e-h scattering in the ambipolar diffusion coeffi- 
cients is also implied in Ref. 8. But unlike the situation 
here, all the Coulomb interaction terms drop out from 
the diffusion coefficient in Ref. 8 completely including 
the coherent part under the ambipolar diffusion approx- 
imation. The absence of the e-h scattering can be eas- 
ily understood, since such scattering represents collisions 
of the electrons and holes, which now are parts of inte- 
gral entities diffusing together under the ambipolar diffu- 
sion approximation. However, the coherent part creates 
a new density- and temperature-dependent environment 
(energy landscape) for the original quasi-particles (with 
unrenormalized energies). They diffuse in the modified 
energy-landscape, resulting in an effective change in dif- 
fusion coefficients, as we will explain in more detail in 
the following. 

To study more quantitatively the effects of the co- 
herent part of the many-body interaction, we examine 
the relative change in diffusion coefficients as defined by 

SD X y = A D X y/D°x Y ‘ In Fi g s * 1-4, plot SDxy (a) 

and D X y (b) with respect to density for all four ambipo- 
lar diffusion coefficients. As model material system, we 
choose GaAs of 8 nm in width. All the material parame- 
ters are standard and will not be listed. Fig. 1 shows the 
familiar density-diffusion coefficient Dnn- I n Fig. 1(b), 
coefficients Dnn (solid lines) and D® NN (dashed lines) 
are plotted versus carrier density at three temperatures. 
The overall feature of the diffusion coefficient is explained 
in detail in a separate paper 9 . We see that the coherent 
many-body effects result in a reduction in the diffusion 
coefficient. The relative change of diffusion coefficient 
is plotted in Fig. 1(a), where we see a diffusion coeffi- 
cient reduction of over 25 percent at 200K. This reduction 
decreases as temperature increases. Similar behavior is 
also observed in Fig. 2 for the mutual-diffusion coefficient 
D T n, which relates carrier density gradient to thermal 
flux Jt- Fig. 3 shows the temperature-diffusion coeffi- 
cient Dtt and D^t ( a ) an d the corresponding relative 
change 8 Dtt- Contrary to the reduction of diffusion co- 
efficients shown in Figs. 1 and 2, we see an increase in 
Dtt , be., Dtt > D^ T . Furthermore, the relative in- 
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FIG. 3: Diffusion coefficient Dtt (b) and its relative change 
5 Dtt (a) versus carrier density at three temperatures as in- 
dicated. Solid and dashed curves in Fig. 3(b) are for Dtt and 
D° r , respectively. 

crease SDtt is much smaller in magnitude than SDnn 
and 8Dtn in Figs. 1 and 2. The change is less than 
10 percent. Similar behavior is shown fn Fig. 4 for the 
mutual-diffusion coefficient Djvt> which describes carrier 
density flux induced by temperature gradient. 

Let us now explain in more detail these figures. We 
begin with the reduction in diffusion coefficients, Dnn 
and D T n, in Figs. 1 and 2, both of which are determined 
by the derivatives with respect to density [see Eqs. (11), 
(12), (17), and (18)]. This reduction can be explained 
by the bandgap renormalization. We note that Dnn de- 
scribes a carrier density flux from high density region 
to low density region. Due to bandgap renormalization 
which increases with carrier density,, the high density re- 
gion has a smaller total bandgap than the low density 
region. This means that a diffusing particle from high 
density region to low density region will have to climb an 
uphill energy landscape due to many-body effects, thus 
leading to a reduction in the effective diffusion coeffi- 
cients. The reduction in the mutual-diffusion coefficient 
Dtn needs slightly different explanation. First, we note 
that Dtn, by definition, describes the thermal flux from 
high density region to low density region. Due to the 
bandgap renormalization, an energy bandgap profile is 
created. A thermal flux is therefore induced from the 
high bandgap (lower density) region to lower bandgap 
(higher density) region to counteract the thermal flux and 
to equilibrate the total energy (bandgap plus the ther- 
mal energy) profile. The situation described in Figs. 3 
and 4 are exactly the reverse of that in Figs. 1 and 
2. The increase in the diffusion coefficients are due to 
the bandgap renormalization that decreases with plasma 
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FIG. 4: Diffusion coefficient Dnt (b) and its relative change 
5 Dnt (a) versus carrier density at three temperatures as in- 
dicated. Solid and dashed curves in Fig. 4(b) are for Dat 
and D% Ty respectively. 

temperature 10 instead of increasing with density as in 
Figs. 1 and 2. The energy landscape reverses from the 
cases of Figs. 1 and 2 and leads to increase of diffusion 
coefficients Dnt an d Dtt • In short, the many-body ef- 
fects on diffusion coefficients lead to a decrease in those 
"diffusTon "coefficTents that are related to the density gra- 
dient [Dnn and Dtn ) a nd an increase in those that are 
related to temperature gradient ( Dtt an d Dnt)- Fi- 
nally, the smaller change in Dtt and Dnt (Figs. 3 and 
4) than in Dnn and Dtn (Figs. 1 and 2) is due to the 
weaker dependence of the bandgap renormalization with 
respect to plasma temperature than to density. 

Another feature of Figs. 1 and 2 is the decrease of 
the relative change SDnn an d SDtn at lower density 
until the carrier density reaches the critical value near 
IxlO 12 cm~ 2 , where electrons become degenerate. This 
decrease is a result of the reduction of the bandgap due 


to bandgap renormalization and to an almost constant 
value of Dnn and d tn • The relative change SDnn 
reaches a minimum around 28 percent for 200K. With 
the further increase of carrier density, SDnn and SD T n 
start to increase as the carriers become strongly degen- 
erate and D° nn and D° TN begin to rise dramatically 9 . At 
high density for lasing over IxlO 12 cm 2 , the diffusion 
reduction is still over 20 percent. The larger values of 
SDnn and SDtn at higher temperature are mainly due 
to the increase of the D and Dj*n with temperature. 

In summary, many-body effects are investigated in a 
semiconductor quantum well where spatial nonunifor- 
mity of densities and temperatures exists along the quan- 
tum well plane. Different roles played by the coherent 
and incoherent parts of Coulomb interaction are ana- 
lyzed. While both coherent and incoherent parts con- 
tribute to the diffusion coefficients of the general two 
component system, the conductivities depend only on 
the scattering part. Even though e-h scattering plays 
an important role in legitimating the ambipolar diffusion 
approximation 7 , we show that the diffusion coefficients 
of the established composite system do not depend on 
the e-h scattering rate. Instead the ambipolar diffusion 
coefficient depends only on the coherent part of the inter- 
action, the bandgap renormalization. We found that the 
coherent many-body interaction leads to a significant re- 
duction of the ambipolar diffusion coefficients Dnn and 
Dtn and an increase in coefficients Dtt and Dnt • 
note that this quite significant change in diffusion coeffi- 
cients, especially in Dnn and Dtn-, should be important 
in describing optoelectronic devices where spatial inho- 
mdgeneities of densities, or plasma temperatures occur. 
Such non-uniformities are quite ubiquitous in high power 
and ultrafast devices, such as lasers and photodetectors. 
Simulation of such devices using the microscopically cal- 
culated diffusion coefficients will be reported elsewhere. 
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